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Downsizing of star-forming galaxies by gravitational processes 
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ABSTRACT 

Context. There is observed a trend that a lower mass galaxy forms stars at a later epoch. This downsizing of star-forming galaxies has been 
attributed to hydrodynamical or radiative feedback processes that regulate star formation. 

Aims. We explain the downsizing by gravitational processes alone, in the bottom-up scenario where galaxies evolve from subgalactic-scale 
objects. 

Methods. Using the theory of Press and Schechter, as a function of galaxy mass, we study the peak epoch of formation of subgalactic-scale 
objects, i.e., gravitational collapse of subgalactic-scale density fluctuation. 

Results. The subgalactic-scale density fluctuation is found to collapse at a later epoch for a lower mass galaxy. The epoch is close to the peak 
epoch of star formation derived from observations. 

Conclusions. The downsizing is inherent in gravitational processes of the bottom-up scenario. Within a region of the initial density field that 
is to evolve into a lower mass galaxy, subgalactic-scale fluctuation is of a smaller amplitude. Gravitational collapse of the subgalactic-scale 
density fluctuation and the subsequent onset of star formation occur at a later epoch for a lower mass galaxy. 
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,^ ' 1. Introduction 



The current paradigm for the formation and evolution of galax- 
ies is the bottom-up scenario based on the cold dark mat- 
ter model. The cold dark matter retains significant small-scale 
fluctuation in the initial density field. Since the first objects are 
formed within dark matter halos that originate in gravitational 
collapse of the initial density fluctuation, these objects are not 
massive. They merge via gravitational clustering of the dark 
matter halos, evolve into low-mass galaxies, and then evolve 
into high-mass galaxies. 

Observationally, there is a trend that is seemingly incon- 
sistent with the bottom-up scenario, i.e., the "downsizing" of 
star-forming galaxies (Cowie et al. [T996 , Heavens et al. 2004 



Thomas et al. 2005; Treu et al. 2005 1. Regardless of environ- 
ment, a lower mass galaxy forms stars at a later epoch. In par- 
ticular, a lower mass galaxy undergoes the peak of star forma- 
tion at a later epoch. 

The existing explanations for the downsizing invoke hy- 
drodynamical or radiative feedback processes that regulate star 
formation (De Lucia et al. 2006 , see also the above references). 
For example, if the galaxy is not massive, supernova explosions 
could expel gas from the galaxy and thereby inhibit star forma- 
tion until the fallback of the gas. However, here we demonstrate 



that the downsizing is in fact explainable by gravitational pro- 
cesses alone, using the theory of Press & Schechter ( 1974) for 
the bottom-up formation of dark matter halos. 

2. Definitions and equations 

The standard values are assumed for the cosmological param- 
eters: the Hubble constant H = 70 km s _1 Mpc -1 (h = 0.7), 
matter density Q. m = 0.3, baryon density Q/, = 0.04, and cos- 
mological constant Qa = 0.7 (Spergel et al. 2006 1. We normal- 
ize the scale factor a(t) to unity at the present epoch. 

The mass density p(x, t) and its average (p) = 3Q.,„Hq/87tG 
lead to the density contrast 5(x, t) = p(x, t)/(p)-l. Here x is the 
position in comoving coordinates and (■) denotes an average. In 
the linear regime, the density contrast grows as 

S(x, t) = D(t)6(x). 

The linear growth factor D(t) is from Heath ( 1977 1: 
Hlda! 
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20 = |><W«" + U - &m - A )/a + Q a] ■ The 
initial density contrast S(x) is related to its Fourier transform 
6(k) as 

1 



S(x) 



(27T)3/2 



6(k) exp(ik ■ x)dk. 
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Since the initial density contrast is random Gaussian, its statis- 
tics are uniquely determined by the initial power spectrum P(k) 
defined as 



P(k) 
A 



(6(k)S(k'y) = (2nY6 D (k - k')P(k). 



(4) 



Here k - \k\ is the wavenumber, 5{k)* is the complex conjugate 
of S(k), and 6u(x) is Dirac's delta function. 

The initial power spectrum originates in the Harrison- 
Zel'dovich spectrum Ck: 



P(k) = T(k) 2 Ck. 



(5) 



We take the transfer function T(k) from Bardeen et al. dl986l l: 
ln(l +2.34 ? ) 



T(k) 



2.34^ 



(6) 



x [l + 3.89^ + (I6.I4) 2 + (5.46 ? ) 3 + (6.1lq) 4 ] 



-1/4 



with q = k/Th Mpc~ l and F = Q m /z exp[-fi ft ( 1 + V2/z/Q,„)] 
(Sugiyama 1995 1. Then, the normalization constant C is set to 
be 7 x 10 4 on the basis of the root-mean-square linear density 
contrast smoothed over a radius Rm = 8 h~ l Mpc observed for 
the present epoch, 0.8 (Spergel et al. 2006 1. 

The smoothing is based on a top-hat window function: 
W M (x) = 3/4nR M for \x\ < R M and for |*| > R M . Here M 
is the mass for both baryonic and dark matter, and R M is the 
corresponding scale: M = AnR^ip) /3. The Fourier transform 
of the smoothed density contrast Sm(x) is 

6 M (k) = W M (k)6(k), (7) 

with W M (k) = 3[sin(kR M ) - kR M cos(kR M )]/(kRM) 3 . This dif- 
fers from the usual Fourier transform of Wm(x) by a factor 

(27T) 3 / 2 . 

3. Basic mechanism 

Consider galaxies that have different masses at an epoch. This 
is consistent with the bottom-up scenario, although high-mass 
galaxies are few if the epoch is too early. The individual galax- 
ies evolve from subgalactic-scale objects. Within such objects, 
star formation begins before the assembly of galaxy mass (see 
De Lucia et al. 2006 1. Using the theory of Press & Schechter 
( 119741 . we discuss that subgalactic-scale objects are formed 
later for a lower mass galaxy. Then, stars are formed later for a 
lower mass galaxy. This is just the downsizing of star-forming 
galaxies. 1 

For convenience, in our discussion, we equate formation of 
an object such as a galaxy with that of its dark matter halo. The 
same convention is used for the mass. 

Press & Schechter 1 9741 considered that, even if objects 
have been formed at a scale, the larger scale density field re- 
mains linear. Regardless of substructure, the object formation 
occurs if the linear density contrast smoothed over that region 
reaches a certain critical value <5 C . Then, at an epoch t, ob- 
jects with mass M exist at all the positions where the smoothed 

1 Neistein et al. 1 2006 1 gave a similar discussion based on a different 
approach. Their preprint was released just before the release of the 
second version of our preprint. 
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Fig. 1. Upper panel: Initial power spectrum P(k). The gray area 
indicates the power removed by the smoothing over the scale 
Rm- Lower panel: Initial density contrast 6(x). The scales Rm 
and Rm' correspond to galaxy masses M and M' (Rm > Rm> for 
M > M'). 



linear density contrast D(t)5M(x) is equal to the critical value 
6 C . The positions with D(t)dM(x) > t> c belong to objects with 
higher masses. Here it is essential to specify the epoch t be- 
cause any object continuously increases its mass in the theory 
of Press & Schechter J1974I . 

The critical density contrast S c is from gravitational col- 
lapse of a homogeneous sphere. In the Einstein-de Sitter uni- 
verse with Q. m = 1 and Qa = 0, the linear density contrast ex- 
trapolated to the epoch of complete collapse is 3(127r) 2 ^ 3 /20 - 
1.69 (Gunn & Gott 179721 This value is used as 8 C . The differ- 
ence from the value for Q m = 0.3 and Q.^ = 0.7 is less than 1% 
up to the present epoch (Percival 2005 1. 

The basic mechanism of the downsizing is illustrated in 
Fig. [2 We compare two regions of the initial density field that 
are to evolve into galaxies with different masses M and M' 
(M > M') at an epoch t. The two regions respectively exist 
at any of the positions where 6m(x) and 5m'{x) are equal to 
6 C /D(t). Superposed on these galactic-scale overdensities that 
have the same height, there are subgalactic-scale density fluc- 
tuations 6<m(x) = 5(x) - 6m(x) and 6<m'(x) = 6(x) - 6m'(x). 
Their variances {6 2 <M ) and {5 2 <M ,) are equal to the powers 
47r J k 2 P(k)dk in the wavenumber ranges k > R^ and k > R^, . 
Since we set Rm > Rm', we have {5 2 <M ) > {5 2 <M ,)- That is, the 
subgalactic-scale fluctuation is of a smaller amplitude for the 
lower mass galaxy. Hence, formation of subgalactic-scale ob- 
jects, i.e., gravitational collapse of the subgalactic-scale fluctu- 
ation, occurs later for the lower mass galaxy. Since the gravi- 
tational collapse induces star formation, the star formation also 
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occurs later for the lower mass galaxy. The lag of star forma- 
tion behind the gravitational collapse is negligible, as far as the 
star formation occurs. This is because the dynamical timescale 
is small in dense dark matter halos that are formed at relatively 
high redshifts. 

4. Analytical calculation 

For formation of subgalactic-scale objects with mass M„ the 
peak epoch f p is calculated as a function of galaxy mass M at 
an epoch t. We set M, = 10 8 or 10 9 Mq as the minimum mass 
for significant star formation. The first stars are formed within 
minihalos with 1O 6 M0, but significant star formation occurs 
after the formation of halos with > 10 8 Mq (Bromm & Larson 
2004] Kitayama & Yoshida l2005t . We also set M = 10 10 , 10 11 , 
10 12 , or 10 13 Mq as the typical galaxy mass. 

The initial density field is smoothed over the scale Rm, that 
corresponds to the mass of subgalactic-scale objects M t . From 
Eqs. @ and @, we have a new power spectrum 



(*) = WuSkTPik). 



(8) 



Eqs. Q, @, @, and yield correlations among the galactic- 
and subgalactic-scale fluctuations 5m(x) and 5<m{x) 



(6 Z M ) = 4tt f eW M (k) z P M XWk, 
Jo 

Jr»oo 
k 2 [\ - W M (k)] 2 PM.(k)dk, 
o 



(9) 



f 

Jo 



I k*W M (W - W M (k)]P M ,(k)dk. 



The correlation (<Jm^<m) has been ignored in Sect.[3] but in fact 
it exists because WM(k) has extended oscillations. Since 6m(x) 
and 5<.m{x) are random Gaussian, their joint probability is 



P(Sm, 6<m) = 
1 

exp 



(10) 



2d 



with d = {S 2 M ){5 2 M ) - {6m$<m) 2 - The probability for 6m(x) is 



p(5 M ) = 



J M'^<M' 
1 



exp 



(11) 



We obtain the probability for the subgalactic-scale fluctuation 
8<m{x) within regions that are to evolve into galaxies with mass 
M at an epoch t. This is the probability under the condition that 
6m(x) is equal to 6 C /D(t): 



PmA<5<m) - 



p[5 M = dJD(t),d <M ] 
p[5 M = SJD(t)] 



(12) 



The conditional probability PmA^<m) is Gaussian. Its average 
and variance are 



(8<m)m,i 



{S 2 M ) D{tY 



{(d<M ~ (5<M)f)M,t = O 2 



(13) 
(14) 
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Fig. 2. Peak redshift z(f p ) of formation of subgalactic-scale ob- 
jects calculated from Eq. dl5l l. The galaxy mass M at redshift 
z(t) is 10 10 , 10", 10 12 , or 10 13 M o . The subgalactic-scale ob- 
ject mass M« is (a) 10 8 Mq or (b) 10 9 Mq. 



The average (5<m)m/ is smaller and hence less important than 
the lcr fluctuation ((r5<M - (5<m)) 2 )^ for the parameter values 
in our calculation. 

The epoch of gravitational collapse of the lcr subgalactic- 
scale fluctuation is used as the peak epoch f p of formation of 
subgalactic-scale objects within regions that are to evolve into 
galaxies with mass M at an epoch t: 



6 C 

D(t) 



+ (8<m)mj + ((<5<m 



<^<M» 2 )^ 



D(t p ) 



(15) 



This equation represents the situation of Fig.^ The lcr fluctu- 
ation corresponds to the typical height for local maxima of the 
subgalactic-scale fluctuation 5<m{x) (Bardeen et al. 1986, Fig. 
2 for r = 0.5). 

The original theory of Press & Schechter i 19741 1 has the 
"cloud-in-cloud" problem. They did not consider correlations 
among different mass scales. In particular, they did not consider 
whether an object belongs to another object with higher mass 
at the same epoch. Judging from the discussion of Nagashima 
( 2001 ), since the correlation (5mS<m) is considered in Eq. dl5> . 
it is rather as accurate as the standard extension of the Press- 
Schechter theory (Bond et al. U991> . There is an exception that 
objects with masses higher than M are not considered in Eq. 
d!5i . but the consequences are not serious as seen below. 

Our results are shown in Fig. [2] where the epochs t and t p 
are converted into the redshifts z(f) and z(f p ). The peak red- 
shift z(t p ) of formation of subgalactic-scale objects is lower for 
lower galaxy mass M. We have thus reproduced the downsiz- 
ing. The dependence of z(f p ) on galaxy mass M is more signif- 
icant for lower galaxy mass M and for higher subgalactic-scale 
object mass M». This is due to the narrower wavenumber range 
Rm ~ k ^ R~m f° r me subgalactic-scale fluctuation. 

The peak redshift z(t p ) of formation of subgalactic-scale ob- 
jects is consistent with the peak redshift of star formation de- 
rived from observations (Heavens et al. 2004 Thomas et al. 
2005 Treu et al. 2005 1, if the subgalactic-scale object mass M, 
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is 10 Mq. The observations have uncertainties and discrepan- 
cies. Our best bet is that galaxies with M 10 12 to 10 13 Mq at 
z(t) = have the star formation peak at z — 4 to 5. 2 These z val- 
ues are close to the corresponding z(t p ) values in Fig.|5Jj. On the 
other hand, galaxies with M % 10 1 1 Mq at z{t) = have the star 
formation peak at z ^ 3. The z values are less than the corre- 
sponding z(t p ) values. This is attributable to the cloud-in-cloud 
problem. The z(t p ) values for Mi 10 11 Mq are overestimated 
because some objects with M S5 10" Mq are substructures of 
galaxies with ^ 1O I2 M0. If we considered this effect, albeit 
necessarily in a complicated manner, the z(f p ) values would be- 
come close to the peak redshifts of star formation. 

The cloud-in-cloud problem also implies that some galax- 
ies belong to clusters of galaxies. Since we ignore this implica- 
tion, we ignore the dependence of the downsizing on environ- 
ment. The dependence is actually negligible (Treu et al. 2005 
see also De Lucia et al. 12006b . Although high-mass galaxies 
with old stellar populations tend to exist in clusters owing to 
the biased galaxy formation (Sect. [5j, high-mass field galaxies 
are equally dominated by old stellar populations. 

5. Discussion 

The downsizing of star-forming galaxies, with a lower mass 
galaxy forming stars at a later epoch, is inherent in gravitational 
processes of the bottom-up scenario. Within a region of the 
initial density field that is to evolve into a lower mass galaxy, 
subgalactic-scale fluctuation is of a smaller amplitude (Fig.[0. 
Gravitational collapse of the subgalactic-scale fluctuation and 
the subsequent onset of star formation occur later for a lower 
mass galaxy. As a function of galaxy mass, we have calculated 
the peak epoch of gravitational collapse of subgalactic-scale 
fluctuation (Fig.|2j. The peak epoch is consistent with the peak 
epoch of star formation derived from observations. 

Our explanation for the downsizing differs from the exist- 
ing explanation for the biasing, where gravitational collapse of 
small-scale density fluctuation occurs earlier if the large-scale 
overdensity in that region has a larger height (Kaiser 1984 
Bardeen et al. 1986 1. We have demonstrated that, if the large- 
scale overdensity has a fixed height, its scale determines the 
epoch of gravitational collapse of small-scale density fluctu- 
ation (Fig.Q. For the biased galaxy formation, the small and 
large scales correspond to galactic and supergalactic scales. For 
the downsizing of star-forming galaxies, the small and large 
scales correspond to subgalactic and galactic scales. The down- 
sizing of star-forming galaxies is thereby independent of the 
biased galaxy formation. 

Thus far, the downsizing has been explained by hydrody- 
namical or radiative feedback processes that regulate star for- 
mation. Although these processes are at work, they would not 
be essential. To the formation and evolution of galaxies, hydro- 

2 Supported by a semi-analytic model on a ,/V-body simulation. The 
star formation peak is at z — 4.5 and 5.0 for ellipticals with M ^ 
10 12 and 10 13 Mq at z(t) = 0, respectively (De Lucia et al. 120061 . 
Such ellipticals dominate the present population of galaxies with M <: 
10 12 Mq. These z values are almost equal to the corresponding z(f p ) 
values in Fig.fJJi. 



dynamical and radiative processes are less essential than grav- 
itational processes that solely explain the downsizing. Also, an 
observation suggests that the downsizing continuously extends 
from galaxies to clusters of galaxies where the star formation 
rate per unit mass at the present epoch is lowest (Feulner et 
al. [2006). Over a wide mass range from galaxies to clusters 
of galaxies, no processes are the same. The exception is grav- 
itational processes, which apply to clusters of galaxies only if 
we consider subgalactic-scale density fluctuation superposed 
on cluster-scale overdensities in the initial field. 

There is observed another trend that a lower mass galaxy 
has a longer history of star formation around its peak (Heavens 
et al. 120041 Thomas et al. 120051 see also De Lucia et al. l2006> . 
This trend, obtained by plotting the star formation rate against 
epoch f, is explainable by the accelerated expansion of the uni- 
verse and the resultant suppression of gravitational clustering 
of dark matter halos. The trend is not significant if the star for- 
mation rate is plotted against the linear growth factor D(t). 

The downsizing is also expected for clusters of galaxies un- 
dergoing galaxy formation. Within a region of the initial den- 
sity field that is to evolve into a lower mass cluster, galactic- 
scale fluctuation is of a smaller amplitude. The galaxy forma- 
tion, i.e., gravitational collapse of the galactic-scale fluctuation, 
occurs at a later epoch. Hence, a lower mass cluster is expected 
to contain younger galaxies. Further studies would be impor- 
tant to understand more about gravitational processes in the 
bottom-up scenario. 
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